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Abstract 



New type of nonlinear (anharmonic) excitations — bushes of vibrational modes — in 
physical systems with point or space symmetry are discussed. All infrared active and 
Raman active bushes for C 60 fulerene are found by means of special group-theoretical 
methods. 



1 Introduction 



Vibrations of many fullerenes and fullerites were investigated by different experimen- 
tal and theoretical methods (see |l| H and references in these papers). Although the 
majority of such studies are based on the harmonic approximation only, some non- 
linear (anharmonic) effects also were discussed in a number of papers. For example, 
combination modes of second order brought about by anharmonicity of interactions in 
C*6o fullerene are discussed and the appropriate lines in infrared transmission spectra are 
reported in ||. The above effects does not exhaust the influence of anharmonicity on 
the fullerene and fullerite vibrational spectra, and we want to consider an application 
for these objects of the consistent group-theoretical approach for studying nonlinear 
vibrations in arbitrary physical systems with discrete (point or space) symmetry de- 
veloped in ^, |5], H [7) . This approach reveals the existence of new nonlinear dynamical 
objects (or new type of anharmonic excitations) in systems with discrete symmetry 
which we call bushes of normal modes. The concept of the bush of modes can be 
explained as follows. 

In the frame of the harmonic approximation, the set of normal modes can be in- 
troduced which are classified by irreducible representations (irreps) of the symmetry 
group G of the considered physical system in equilibrium. In this harmonic approxi- 
mation, normal modes are independent of each other, while interactions between them 
appear when some anharmonic terms in the Hamiltonian are taken into account. Let 
us note that a very specific pattern of atomic displacements corresponds to each nor- 
mal mode. As a consequence, we can ascribe to a given mode a definite symmetry 
group Gd which is a subgroup of the symmetry group G. The group Gd is a symmetry 
group of the instantaneous configuration of our system in its vibrational state. 

Let us excite at the initial instant to only one, arbitrarily chosen mode which will 
be called the root mode. We suppose that all other modes at the initial moment have 
zero amplitude. Let the symmetry group Gd and irrep r correspond to this root 
mode. Then we can pose the following question: to which other modes can this initial 
excitation spread from the root mode? We will refer to these initially "sleeping" modes, 



belonging to the different irreps (j 7^ 0), as secondary modes. 

A very simple answer to the above question was found in 0, [|]. It turns out that 
initial excitation can spread from the root mode only to those modes whose symmetry 
is higher than or equal to the symmetry group Go of the root mode. We call the 
complete collection containing the root mode and all secondary modes corresponded 
to it a bush of modes. Since no other modes are excited, the full energy is trapped 
in the given bush. As a consequence of the above idea, we can ascribe the symmetry 
group Gd (remember that this is a group of the root mode) to the whole bush, and in 
this sense we can consider the bush as a geometrical object. 

It was proved in @, [5], || that all modes belonging to a given bush B [Go] are 
coupled by force interactions. It is very important that the structure of a given bush 
is independent of the type of interactions between particles of our physical system. 

A bush of normal modes can be considered as a dynamical object, as well. Indeed, the 
set of modes corresponding to a given bush B [Go] does not change in time, while the 
amplitudes of these modes do change. We can write exact dynamical equations for the 
amplitudes of the modes contained in the bush B [Go], if interactions between particles 
of our physical system are known. Thus, the bush B [Go] represents a dynamical system 
whose dimension can be essentially less than that of the original physical system. 

The above properties of bushes of normal modes can be summarized in the following 
manner. A normal mode represents a specific dynamical regime in a linear physical 
system which, upon being exciting at the initial instant to, continues to exist for any 
time t > to. Similarly, a bush of normal modes represents a specific dynamical regime 
in a nonlinear system which can exist as a certain object for any time t > to- 

2 Some mathematical aspects of bushes of normal 
modes 

Let us examine a nonlinear mechanical system of N mass points (atoms) whose Hamil- 
tonian is described by a point or space group G. Let three-dimensional vectors Xj(£) (i = 



1,2,..., N) determine the displacement of the i-th atom from its equilibrium position 
at time t. The 3 x iV-dimensional vector X(£) = {xi(t), x 2 (t), . . . , xjv(t)}, describ- 
ing the full set of atomic displacements, can be decomposed into the basis vectors 
(symmetry-adapted coordinates) of all irreps 1} of the group G contained in the me- 
chanical represent ationQ 

xw = E^rf = E^, (i) 

P j 

Here iff is the z-th basis vector of the n., -dimensional irrep Fj. The time dependence 
of X(£) is contained only in the coefficients Hji{t) while the basis vectors are time 
independent. Thus, a given nonlinear dynamical regime of the mechanical system 
described by the concrete vector X(t) can be written as a sum of the contributions Zlj 
from the individual irreps 1} of the group G. 

Each vibrational regime X(t), can be associated with a definite subgroup Gd {Go C 
G) which describes the symmetry of the instantaneous configuration of this system. 
Now the following essential idea is proposed. The subgroup Gd is preserved in time; 
its elements cannot disappear during time evolution except in the case of spontaneous 
breaking of symmetry which we will not consider in the present paper. This is the 
direct consequence of the principle of determinism in classical mechanics. 

Introducing operators g G G acting on the 3iV-dimensional vectors X(t), which 
correspond to elements g G G acting on three-dimensional vectors Xj(£), we can write 
the above condition of preservation of Gd as a condition of invariance of the vector X(t) 
under the action of the elements of the group Gd' 

g X(t) = X(t), g G Gd- (2) 

Combining Eqs. (|l|) and (|2|) one obtains (for the details see 0) the following invariance 
conditions for individual irreps Ff 

(Fj | Gd) Cj = Cj. (3) 

1 Considering vibrational regimes only, we can treat F as a 3iV — 6 vibrational representation of the 
group G. 



Here 1} j Gd is the restriction of the irrep 1} of the group G to the subgroup Gp, 
i.e. the set of matrices of /} which correspond to the elements g G Gd only. The rij- 
dimensional vector Cj in Eq. (|3]) is £/ie invariant vector in the carrier space of the irrep 7} 
corresponding to the given subgroup Gd C G. Note that each invariant vector of a 
given irrep 7} determines a certain subspace of the carrier space of this representation, 
and the total number of arbitrary constants upon which the vector depends is equal 
to the dimension of this subspace. If in solving Eq. (|3|) we find that Cj ^ 0, then the 
irrep 7} does contribute to the dynamical regime X(£) with the symmetry group Gd- 
Moreover, the invariant vector Cj determines the explicit form of the mode of the 
irrep 7} belonging to the bush of modes associated with the given nonlinear dynamical 
regime. 

We will illustrate the general statements of bush theory with the Cqq fullerene with 
buckyball structure and with the icosahedral symmetry group G = Ih- There are 
10 irreducible representations of dimensions 1 (A g , A u ), 3 (Fi g , F\ u , F 2g , i*2 U ), 4 (G s , 
G u ) and 5 (H g , H u ) associated with the group Ih- The infrared (IR) active modes 
belong to the irrep F\ u , and modes which are active in Raman (R) experiments belong 
to irreps A g or H g . We found all bushes of modes for Cqo fullerene. There are 22 
different bushes for this fullerene. Let us consider the bush B7 corresponding to the 
symmetry group Gd = C 5v C Ih- Only four irreps A g ,H g ,Fx u and F 2u contribute 
to it (the appropriate invariant vectors are zero for all other irreps of the icosahedral 
group G = I h ) Q: 

B7 : [symmetry C 5v ] : A g (a) - I h , H g (a, 0.577a, 0, 0.516a, -0.258a) - D 5d 

(4) 

F lu (0,0,a)-G 5 „ F 2u (a, 0.258a, 0.197a) -C 5v . 

The arbitrary constants entering into the description of different invariant vectors are 
not connected with each other. As all invariant vectors listed in Eq. (Q) are one- 
parametric (their arbitrary constants are denoted by the same symbol a only for clar- 
ity), it is clear that the bush B7 depends on four arbitrary constants (one constant 
2 Since some elements of the matrices of multidimensional irreps of the group G = Ih are irrational 
numbers, we keep only three digits after the decimal point when we write the invariant vectors. 



for each of the four irreps). The structure Eq. (H) of the bush B7 shows that there 
exist only four contributions Aj to the appropriate dynamical regime X(£). We denote 
themQ as 2lL4 g ], A[H g ], A[Fi u ], A[F 2u ]. The invariant vectors listed in Eq. ([|) permit 
us immediately to write the explicit form of the dynamical regime X(£) correspond- 
ing to the bush B7 by replacing the arbitrary constants with the four functions of 
time /i(t) , v(t) , j(t) and £(t): 

X(t) = A[A g ] + A[H g ) + A[F lu ] + A[F 2u ] = 

ti(jL)<p[A g ] + v{t){<p x [H g \ + 0.577cp 2 [H g ] + 0.516<p 4 [# s ] - 0.258(^ 5 [H g ]}+ (5) 
j{t)(p 3 [F lu ] + e(t){^[F 2 J + 0.258ip 2 [F 2u ] + 0. 197(^3 [F 2u ]}. 

Eq. (H) is a consequence of the relation of the group G and its subgroup Gd only, and 
now we must take into account the concrete structure of our physical system to find 
the explicit form of the basis vectors ip\ of the irreps entering into Eq. (||). They can 
be obtained by conventional group-theoretical methods, for example, by the projection 
operation method. The basis vectors of the irreps determine the specific patterns of 
the displacements of all 60 atoms of the 6*60 fullerene structure. 

It is important to note that each of the irreps A g ,H g ,F\ u and F 2u is contained 
in the vibrational representation of the 6*60 fullerene several times, namely, 2, 8, 4 
and 5 times, respectively. (These numbers are equal to the numbers of fundamental 
frequencies of normal modes associated with the considered irreps). As a consequence, 
we must treat the time-dependent coefficients in Eq. (||) as vectors of the appropriate 
dimensions. Because of this we ascribe a new index (k) to the basis vectors determining 
the number of times (rrij) the irrep 1} enters into the vibrational representation. Each 
contribution Aj "splits" into rrij copies Ajk, where k — 1, 2, . . . , rrij and, therefore, 

x(i) = EA = HEA*). (6) 

j j k=l 

For the case of bush B7 we have A[A g ] = A x [A g ] + A 2 [A g ], A[F lu ] = Ai[F lu ] + 
A 2 [F lu ] + A 3 [F lu ] + A A [F lu ], etc. 



3 Hereafter we write the symbol j of the irrep generating the contribution Aj in square brackets 
next to symbol 



The bush B7 in the Cqq fullerene structure forms a 19- dimensional dynamical ob- 
ject: its evolution is described by the dynamical variables listed below as compo- 
nents of the four vectorial variables /x(t) , v(t), "f(t) and £(£): fj,(t) = [[ii(t), ^(t)], 
u(t) = [u 1 (t),..., u 8 (t)), 7 (t) = [ Tl (t), . . . , T4 (t)], = [&(*), • • • , &(*)]• 

Thus, although only four of the ten irreps contribute to the bush B7, its dimension 
is equal to 19 because several copies of each of these four irreps are contained in the full 
vibrational representation of Cqq fullerene. We cannot predict the concrete evolution 
of the amplitudes of the bush modes without specific information of the nonlinear 
interactions in the considered physical systems, but we can assert that there does exist 
an exact nonlinear regime which involves only the modes belonging to a given bush. 

3 Optical bushes for Cqq fullerene 

As was already noted, there are 22 bushes of vibrational modes for Cqq fullerene struc- 
ture. Five of them are infrared active and six are Raman active. We call these bushes 
by the term "optical". The root modes of the optical bushes belong to the infrared 
active irrep F lu or to the Raman active irreps A g and H g . We want to emphasize that 
some modes associated with the irreps which are not active in optics can be contained 
in a given optical bush. 

All optical bushes with their symmetry groups (in square brackets), numbers of 
irreps contributing to them and their dimensions (in parentheses) are listed below. 

Infrared active bushes: 
B7 [C 5v ] (4, 19); Bll [C 3v ] (6, 31); B15 [C 2v ] (7, 46); 

B19 [C a ] (9, 89); B22 [d] (10, 174). 
Raman active bushes: 
Bl [4] (1, 2); B4 [D M ] (2, 10); B5 [D u ] (3, 16); 

BIO [D 2h ] (3, 24); B16 [C 2h ] (5, 45); B20 [Cj (5, 87). 

Supposing that nonlinearity of the considered system is u>ea£j] we can estimate the 
relative values of the contributions from different irreps to a given bush. For example, 
4 According to results obtained in the paper || this hypothesis valid for Ceo fullerene vibrations. 



for above discussed infrared active bush B7 we have 

A[F lu ](root) = 0(e); A[F lu ](secondary) = 0(e 3 ); 
A[F 2u ] = 0(e 3 ); A[A g ] = 0(e 2 ); A[H g ] = 0(e 2 ). 

Here e is an appropriate small parameter characterizing the value of the root mode. 

Thus, in the case of weak nonlinearity, the contributions of different irreps can be of 
essentially different value. This property seems to be important for the interpretation 
of the vibrational spectra of bushes of modes. 

4 Conclusion 

In the present paper, we consider a new type of possible nonlinear excitations — bushes 
of normal modes — in vibrational spectra of fullerenes and fullerites, using as an ex- 
ample the 6*60 buckyball structure. We believe that special experiments for revealing 
the bushes of vibrational modes in their pure form will be important for further elu- 
cidation of the role of these fundamental dynamical objects in various phenomena in 
fullerenes and fullerites. It seems that such experiments may be similar to those by 
Martin and others reported in || . However, unlike these experiments, we must use the 
monochromatic incident light with frequency close to that of the root mode and with 
polarization along the symmetry axis of the chosen bush. 

The first-principle calculations are desirable for obtaining the coefficients of the 
anharmonic terms in the Cqq fullerene for a more detailed description of the bush 
dynamics. 

The concept of bushes of normal modes and the appropriate mathematical methods 
for their analysis are valid for both molecular and crystal structures. Such a possibility 
can simplify the assignment of the different optical lines in fullerites brought about by 
both intra- and inter-vibrations of the Cq$ molecular clusters. 

It will be very interesting to study interactions between bushes of vibrational modes 
and electron subsystems in fullerenes and fullerites. 
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